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1. INTRODUCTION 



The Calogero-Sutherland model has recently attracted some attention mainly because, 
in spite of its simplicity, it yields nontrivial results which contribute to shape our under- 
standing of fractional statistics Its most remarkable property is that its ground state 
is given by a Jastrow wave-function which is the one dimensional analogue of the Laughlin 
wave- function. The excitations are described by quasiparticles which carry a fraction of the 
quantum numbers of the fundamental particles. The wave functions are simple enough to 
allow the computation of physical quantities such as the dynamical correlation functions. 
Recently, it has become clear that the spin version of this model is closely related to the 
known spin chains such as the XXX chain. The wave functions of the spin models can 
be obtained as simply as those of the Calogero-Sutherland model by diagonalising simple 
differential operators known as the Dunkl operators 0. 

In these lectures, I restrict to the simplest models and I introduce some techniques use- 
ful in their study. These lectures are based on work done in collaboration with D.Bernard, 
M.Gaudin and D.Haldane || and with F.Lesage and D.Serban J|]. 

In the first part, I review Sutherlands method to diagonalise the Calogero-Sutherland 
Hamiltonian ||. I then give the expression of the form factor of the density operator. On 
this example, I describe the fractional character of the quasiparticles which propagate in 
the intermediate states. 

In the second part, I review the spin generalisation of the Calogero-Sutherland model 
and I diagonalise their Hamiltonian using the Dunkl operators. I then exhibit a repre- 
sentation of the Yangian which commutes with the Hamiltonian. This is achieved by 
"quantizing" the spectral parameters of a monodromy matrix obeying the Yang-Baxter 
equation. The consistancy of this procedure requires that the spectral parameters (the 
Dunkl operators) obey the defining relations of a affine-Hecke-algebra. Finally, I obtain 
a representation of this algebra which degenerates to the Dunkl operators using some 
operators defined by Yang in his study of the ^-interacting gas . 

In the last part, I review the long range interacting spin model known as the Haldane- 
Shastry chain Although closely related to the Calogero-Sutherland model, this chain 

is more difficult to study because the Dunkl operators cannot be used to obtain the wave- 
functions in a straightforward way. I exhibit a representation of the affine Hecke algebra 
in terms of parameters. This representation becomes reducible at the special point where 
it commutes with the Haldane-Shastry Hamiltonian. I then use a correspondance with the 
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Calogero-Sutherland models to obtain the eigenvalues and some of the eigenvectors of this 
Hamiltonian. 

The most unphysical feature of the models discussed here is the long range character 
of the particle particle interaction which behaves as 1/x 2 in the thermodynamical limit. 
More realistic models with short range potential (1/ sinh(x) 2 ) can be defined || but are 
more difficult to study (see [[H| for some recent progress). Some aspects of the conformal 
limit of the Haldane-Shastry chain are also discussed in [Jlljj . 



2. THE CALOGERO-SUTHERLAND HAMILTONIAN 

This model describes particles on a circle interacting with a long range potential [3,4]. 
The positions of the particles are denoted by Xj, 1 < i < N, < xi < L and the total 
momentum and Hamiltonian which give their dynamics are respectively given by: 

N 1 d 



< * n nnr ■ 
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From now on, we shall work on the unit circle and set 9j = 27r^- and Zj = e l9j . The wave 
functions solution of the equation Hip = Eift can be given the following structure: 

m = (p(9)A ,3 (9) (2.2) 

with 

A(9) = ]J S m( 6 -^A (2.3) 

i<3 ^ ' 

$ is a symmetric polynomial in the variables Zj = e ldj and zj 1 . 

To understand the simplicity of the spectrum, one must write the effective Hamiltonian 
H = A-PHAP acting on 0. 
Following Sutherland 0, we obtain: 

The remarkable property of this Hamiltonian is that acting upon symmetric polynomials 
of a given homogeneity in the variables Zj it is realised as a triangular matrix. So, it's 
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eigenvalues can be read on the diagonal of the matrix and there are simple algorithms to 
find the eigenvectors. 

More precisely, let us define the following basis of symmetric polynomials indexed by 
a partition: Ai > A2... > A^v > 0. 

which is the sum over distinct permutations of the monomial z^ z 2 2 ...z N N . 

First, it is easy to see that the subspace of polynomials of a given homogeneity (| A | = 
Y2iLi ^i) i s preserved. Then, inside this subspace, we can define an order on the partition 
by saying that A > p, if Ai > /zi, Ai + A2 > pi + /U 2 , Ai + ... + Xn > A*i + ^2 + ••• + A* at. 
It follows from Sutherland's argument that: 

Hm\ = ^2 c ^ m v ( 2 - 6 ) 
that is to say H is a triangular matrix. The eigenvalues are given by the diagonal elements: 

JV 



i=i 



E x = C XX = Y, (\ + (N- - ((AT - z)/3) 2 (2.7) 
The momentum is given by: 

Px = J2 ^ + ( N ~ 0)9) (2-8) 



N 



i=l 



This set of states is not complete, but the complete set can be obtained by multiplying 
these wave functions by (ziZ2--.zn) p where p is a positive or negative integer. Taking this 
into account, a complete basis is given in terms of N momenta ki = Xi + (N — i)(3 + p and 
the energy and momentum are given by: 



P\ki...k N ) = 

' N \ 

J2 k i )\ki-k N ) 

yi=l J 



(2.9) 

H |fci...fcjv) 



This looks like a free particle spectrum, in particular for (3 = 0,1, the Hamiltonian is a 
pure kinetic energy term and it reduces to the usual boson and fermion description of the 



states. To understand the difference, let us consider the case of integer values of f3. One 
sees from the expression of the k{ that they are integers and obey the constraint: 

k-k i+1 >p (2.10) 

If (3 = 1, this is just the Pauli-exclusion principle, but if (3 = 2, 3, ... the Pauli principle is 
replaced by a stronger exclusion principle. 



2.1. Particle hole excitations 

We shall now describe the excitations of the ground state of such a system and show 
that they can naturally be described in terms of quasiparticles. To illustrate this, we shall 
consider the matrix elements (form factor): 

(a | p{x) | 0) (2.11) 

of the density operator 

N 

p(x) = J2H0i-x) (2.12) 

i=l 

where the state | > stands for the ground state and | a > for some excited state. The 
results which we present here were motivated by the evaluation of the density density 
correlation function obtained by Simons, Lee and Altshuler in the matrix models |12|] (see 
also §). 

The question we address is for which states | a > is the matrix element non zero? 

In the case f3 = 1 which corresponds to free fermions, it is easy to answer using second- 
quantization arguments. The operator p{x) can be represented as i(j + (x)ip(x) where i^(x) 
is a fermion field. Acting on the vacuum, the operator ij){x) anihilates a particle from the 
Fermi sea (| k \ < kp where kp is the Fermi momentum) and the operator i(j + (x) creates a 
particle out of the Fermi sea. The states | a > are therefore described by (A^ — 1) momenta 
ki less than the Fermi momentum kp and 1 momentum k^ larger than kp or equivalently 
in terms of a hole of momentum less than kp and a particle of momentum larger than kp. 

For f3 arbitrary integer, the calculation is based on the properties of the wave functions 
of H (the Jack polynomials) established by Macdonald and Stanley |T3| |TJ[ . The final result 



is however simple enough to be described here. Let us view the ground state as a Fermi-sea 
filled with N particles subject to the constraint \ki — kj\ > (3 : 

k i = (3(i-^—\ l<i<N (2.13) 



The Fermi momentum is given by the largest value of ki : 



N - 1 

k F =(3^—. (2.14) 



This is represented on Figure 1 a) for (3 = 2, N = 5. The ones stand for the integers ki 
occupied by a particle and the for the remaining integers. 

i i 
i i 

a) 00 010101010 10000 

i i 

b) 1 X 1 X 1 1 1 



F F 



Fig. 1: The Fermi sea 



It turns out that the states | a > which propagate are those for which (N — l) momenta 
are inside the Fermi-sea \ki\ < kp and one momentum is out k^ > kp (see Figure 1 b) 
). Such states can be completely described by (3 holes of momentum less than kp and a 
particle of momentum larger than kp. A hole corresponds to a sequence of (3 consecutive 
integers not occupied by a momentum (represented by a cross on the Figure). 

The outcome of the computation [Q (These results have also been conjectured in 
15|j and |16[]) is that these holes are eigenstates of a Calogero-Sutherland Hamiltonian 
with coupling constant (3~ x and have a mass equal to 4 the mass of the particles. In the 
thermodynamic limit the form factor is equal to : 

(a | p(0) | 0) = (w* - 1) ^ f[(l - vf)^ J] \vi - vj\^ ^; Vi - pW \ (2.15) 



where v i and w are respectively the rapidities of the holes and of the particle describing the 
intermediate state. Other form factors such that (A | ij)(x) \ 0) confirm this description of 
the excitations. Although the description given here is reminiscent of a second quantized 
formalism, I am not aware that such a formalism exists. 



3. THE SPIN CALOGERO-SUTHERLAND MODELS 

A step towards the unification of the Calogero-Sutherland models and the integrable 
spin chains consists in introducing spin generalizations of the model studied in the last 
section. 

The Hamiltonian we consider is 

N , o \ 2 

j=l V OZ lS iftj ( Z i ~ z j) 

The particles have an internal spin degree of freedom and the permutation Pij exchanges 
the spins of the particles i and j. 

The eigenstates have the following structure: 

iP (z u a,) = $ ( Zi , a t ) \\ ( Zi - Zj f (3.2) 

i<j 

where the wave function (j> crj) is completely antisymmetric under the simultaneous 
permutations of the spins and coordinates. 

For our purpose, it is more convenient to look at the effective Hamiltonian acting on 

It turns out that this Hamiltonian can be diagonalized in a similar way as the one intro- 
duced in the last section. It has the same additive spectrum but the selection rules for the 
momenta differ from (|2.10| ). 

To do this, let us define the following Dunkl operators: 

* = - & -^Err (*« " i) +^7^7 - 1) (3-4) 

where permutes the coordinates z% and 2j : 

[^,z fe ]=0 for k^i,j , . 

K-z=zK [ } 

Some motivation for the form of these operators and the algebra which they obey will 



come later in the text (see also |L7j [[18]] ) , but first, let us see how they can be used to 
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solve the model. We shall show that the diagonalization of H can be reduced to the 
simultaneous diagonalization of the Dunkl operators di. Moreover in a convenient basis 
the di are represented by triangular matrices. 

The di obey the defining relations of a degenerate affine-Hecke algebra: 

a) [di, dj] = 

b) [K ii+1 ,dk}=0 if k^i,i + l (3.6) 

c) K iji+1 di - d i+1 K iti+1 + (3 = 

From these relations, one deduce that the quantity: 

N N 

D(u) = J] (u - di) = u N + c k u N ~ k (3.7) 

i=l k=l 

commutes with the permutations K^. 

By computing explicitely dj = cf — 2c2, one sees that it is equal to the Hamiltonian 
H except that the P^- are replaced by —Kij. Since D(u) commutes with the permutations, 
we can restrict to the space of antisymmetric wave functions <fi(zi,(Ji) . In this subspace 
— K^ and Py coincide and dj is equal to H. Thus, instead of H, it is easier to simul- 
tanesouly diagonalize the di 's. 

If we consider the explicit expression ( |3.4| ) of the d^'s, one sees that they preserve the 
space of homogeneous polynomials of a given degree in the variables z%. It then follows 
from simple arguments that these operators are realized by triangular matrices in the 
monomial basis Their eigenvalues are equal (up to a permutation) to 

the components of the multiplet: 

5 l = \ l + f3{N -i) (3.8) 

The eigenvectors of dj corresponding to a given eigenvalue form a representation of the 
permutation algebra. The wave functions (zi, <7j) are obtained by antisymmetrizing these 
eigenvectors with an arbitrary spin function X (a^) . In particular, if we antisymmetrize 
the wave functions, we recover the Jack polynomials of the last section (at coupling + 1) 
multiplied by the Vandermonde determinant Yli<j ( z i ~~ z j) ■ The mathematical properties 
of the Jack polynomials can probably be extanded to these wave-functions. 
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3.1. Dunkl operators and monodromy matrices 

It is clear from the explicit construction of the wave functions of H that it's spectrum 
is highly degenerate and we shall now interpret these degeneracies from a symmetry prin- 
ciple. For this, we need to introduce the so-called monodromy matrix T(v) solution of the 
equation: 

R ab T a T b = T b T a R ab (gg) 

where T a in the operator valued matrix T (u a ) <8> 1, T b is 1 <E> T (ub) and R ab is given by 

R ab = u a - u h + (3P ab (3.10) 

where P ah permutes the two components of the tensor product V a <8> V b . 

Consistency of the above relations requires that R ab is a solution of the Yang-Baxter 
equation in V a (g> V b £g> V c : 

R ab R ac R bc = R bc R ac R ab 

Assume that T(u) has a series expansion in 1/u, the coefficients of this expansion define 



an algebra which is called the Yangian |19[]. In this section, we exhibit a representation 
of the Yangian which commutes with H. For this, let us define the so-called "L operators 
" defined as : 

L l = u-d l -pP at (3.12) 

where P ai permutes the components of the tensor product V a <8> V 1 and di is a coefficient 
which we identify with the Dunkl operator (|3.4|) . This is consistent since the d^s commute 



among themselves and with the permutations Pij. A well known representation of the 
monodromy matrix T(u) is then given by ||20|| : 

T a (u) = L 1 (u)L 2 (u)...L N (u) (3.13) 

Now, let us act with this monodromy matrix on the space of wave functions (j)(zi,<Ti) . 
The coordinates Zi are permuted by the Kij and the spins o~i are permuted by the P^. 
We shall see that the antisymmetric wave functions (f>(zi,ai) are preserved by the action 
of the monodromy matrix T{u). To show this, we define a projection IT which substitutes 
—Pij for to the right of an expression. 

The subspace of antisymmetric wave functions is preserved if the following condition 
is satisfied: 

U((K lJ +P lJ )T(u)) = (3.14) 
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which is equivalent to the condition: 



n {{K ii+1 + P u+1 ) U{u)L i+l {u)) = (3.15) 

when we expand this relation in powers of u, the coefficients of u and u 2 give relations 
which are trivially satisfied. Let us here consider the coefficient of 1. We omit the IT in 
front of the expression: 

= (K 12 + P 12 ) (<M 2 - Pd,P a2 - f3d 2 P al + f3 2 P al P a2 ) (3.16) 

using ( |3.6| ) and the fact that [Ki 2 , did 2 ] = we obtain: 

—(3PalPa2Pl2 + fiP\ 2 P a \P a2 

-d 2 P 12 P a x-d x P X2 P a2 (3.17) 
+P a i {d 1 P l2 -(3)+ P a2 (d 2 P 12 +P) = 

which is easily shown. Notice that the unexpected term equal to (3 in ( |3.6| )c) is necessary 
for the proof of this relation. 

It follows from this analysis that the monodromy matrix T{u) generates a representa- 
tion of the Yangian algebra in the subspace of antisymmetric wave functions (z{, <7j) . It is 
also easy to see that this representation commutes with the Hamiltonian H and therefore, 
the Yangian is the symmetry algebra which explains the degeneracies of H. 

3.2. Yang's representation of affine-Hecke algebras 

In this section, we establish a relation between the Dunkl operators ( |3.4| ) and some 
operators defined by Yang in his study of the S interacting gas (see also [^TJ for a discussion 
in the elliptic case). 

To make this relation more transparent, we shall consider the case of the affine-Hecke 
algebra defined by : 

«) [yi:Vj] = o 

b) = if + 1 (3.18) 

c) giVi - yi+ig' 1 = o 

and the gi obey the relations of a Hecke algebra : 

[9i,9j] =0 if \i-j\>2 

9i9i+i9i = 9i+i9i9i+i (3.19) 
(9i - Q) [9i + q' 1 ) = o 
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This algebra replaces the permutations and the d^s if one repeats the arguments of the last 
section replacing the rational solution (|3.10| ) of the Yang-Baxter equation ( |3.11| ) with a 



trigonometric solution. This is considered in H . This algebra is also studied in the works 



of I.Cherednick ||17|| . Here we shall content ourselves to obtain a natural representation of 
the affme-Hecke algebra ( |3.18| ) which in some limit degenerates to the Dunkl operators 

For this, let us consider some realization of the Hecke-algebra (|3.19|) acting in the 
tensor product V® N where the operator gi acts in Vi (g> V*+i. To ease the forthcoming 



discussion, we rename gi : gi^+i- Following Yang's argument ||22|| . we shall exhibit a set of 
operators which commute together. Then we shall see that in some limit these operators 
can be identified with the i/j. 

The starting point is the trigonometric form of the L operator |2(J given by: 



Vi 



LUv)=[-g ia -g- 1 )K at (3.20) 



where g a i acts in V a ® Vi and K ai permutes the two components of V a ® Vi. Using this 
form we can construct a monodromy matrix solution of ( |3.9| ) given by : 



T a {u) = L 1 {u)L 2 (u)...L N {u)S a (3.21) 

where S a is a matrix acting in V a and such that [R ab , S a ® S ] = 0. It then follows from 
the Yang-Baxter equation ( |3.9| ) that the trace of this monodromy matrix T{y ) defines a 
commuting set of transfer matrices: 

[f(v),f(w)] = Vu, (3.22) 

To define the trace of the monodromy matrix, we consider T a as a matrix acting in 
V a and take its trace in the usual sense. 

Hence, the operators jji = T (v^ form a commuting set . These operators play an 
important role in the study of the S interacting gas [2^] . They can be explicitely computed 
using the fact that: 

Li {Vi) = {q- q- 1 ) K ai (3.23) 

If one then writes T (v^ in the form T (v^ = tr Ki a 6 where 9 is an operator that contains 
only the permutations Kij with i,j^a and uses the following property of the trace: 

tr K ai 6 = 6 (3.24) 
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one finally obtains : 

Vi = Xii + iXa + 2...XnsrS l Xii...Xii-i (3.25) 
where the operators Xij are given by: 

Xa= {^-ga-g-^K^ (3.26) 

The operators jji constructed in this way obey the relation Q3.18|) a) but not ( |3.18|) b) c). It is 
then not difficult to see that these relations are satisfied in the limit <C v± <C v<i--- <C vn- 
In this limit the operators yj are given by 0] [23] : 



Vi = g ii l 1 K ii+1 ...g il ^K iN S l K li g li ...K i _ 1)i g i _ lji (3.27) 

To make contact with the Dunkl operators di, we must use the following representation 
of the Hecke algebra acting in the space of polynomials in N variables Zi : 

g ii+1 = qK ii+1 - (q - q' 1 ) (K ii+1 - 1) (3.28) 

Zi — Zi+i 

and S l is defined by: 

(5V)(z ls z N ) = f(z l: ...,tzi, ...,z N ) (3.29) 
If one sets q = t@ and let t — > 1, the operators yi tend to the di defined in (|3.4|) . 



4. THE HALDANE SHASTRY CHAIN 

Among this class of models, the Haldane Shastry chain pj] M is the closest to the known 
integrable spin chains such as the XXX chain ||22|| . It is also the less well understood. Here, 
we shall present these chains from the point of view of the Dunkl operators. 

The Hamiltonian is the (3 = oo limit of (|3.1|), it is equal to : 



H = Y—^1 (4.1) 

^sin 2 ^-^ 



where the variables 9j take the values < j < N. The main difficulty comes from the 
fact that this Hamiltonian, unlike ( |3.1|) , cannot be obtained from the static limit of the 
Dunkl operators as ^2 because in this limit this quantity is a number. In this section, we 
define a representation of the degenerate affine Hecke algebra (|3.6| ) labelled by N complex 
numbers u>k- When the ujk are arbitrary, the representation is irreducible, when they are 
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equal to exp(i^^), it is reducible and one can find the Haldane Shastry Hamiltonian in 
its center. We shall then use a correspondence with the /3 = 1 representation ( |3.4| ) to 
diagonalise the Hamiltonian (|4.1[) . 

Consider the representation of the degenerate affine Hecke algebra ( p.6| ) given by: 

It acts on the Hilbert space defined by the permutations : 

H N = {\Px...P N > /Pi e Z*, P ^ Pj for all i, j} (4.3) 
in the following way : 



K^j | . . . Pi . . . Pj . . . ^> | . . . Pj . . . P^ ... ^> 



(4.4) 



z j \P 1 ...P N > = Up.\P X ...P N > 

In the case where the ui^s are arbitrary, there is no reason for this representation to be 
reducible. On the other hand, if Uk = e l ^r ? the operator C defined by: 

C\P j > = \(P j + l)> (4.5) 

obviously commutes with the dj's and the Kij. So, in this case the representation is 
reducible and one can verify that the Hamiltonian given by : 

i<j j 

is also in the center. No systematic way to obtain the operators which commute with the 
K^a and the d^s is known. In what follows we restrict to the case where Uk = exp(i^£). 



The Hamiltonian (|4.6| ) can easilly be related to the Haldane Shastry Hamiltonian and 



we shall now indicate how to obtain its eigevalues and some of its eigenvectors. We shall 
construct a class of states on which the two representations of the di ( |3.4j ) (for (3 = 1) and 
(fO|) act in the same way. On these states, any conserved quantity of the spin chain has 
its analogue in the dynamical model and can thus be diagonalised . 
For M < N let us define the Hilbert space 1~Lm 

Km = {\Px...Pm > /Pi e Z N , Pi / Pj for all Lj} (4.7) 



12 



which is uniquely embedded in Hn if one requires that 



K ii \P x ...P M >= \P 1 ...P M > fori,j>M (4.8) 
In Hm we consider the following states : 

\a>= V(uj p \...,uj p ™)\P 1 ,...,P m > (4.9) 

where : 

Ufa) = QMjlizi - Zj ) (4.10) 

i<3 

and ty(zi) is a polynomial of degree between 1 and iV — 1 in the variables Z{. On these 
states, di = di — for 1 < % < M. It results from the fact that a polynomial P(z) of 
degree between 1 and N — 1 satisfies the relation : 



(zd z - ^±l)P(z)\ x=Uh = -^T P (^) ( 4 - n ) 
2 v k Uk ~ Ul 



The same relation enables to show that in this subspace H acts as Y^k=\{^k + 1/2 + 
N/2)(d k + 1/2 - N/2) and to find its spectrum: 

M 

E = J2 <^ with e ( fk ) = ( fk + X /2) 2 - N 2 /4 (4.12) 
k=i 

fk + ^-f - are the degrees of the highest monomial of *&(zi). Therefore, the fk satisfy the 
inequalities: 

N - 1 N - 1 

7T- <fl< /2- < - 1 (4.13) 
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